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Is hydrodynamics relevant to RHIC collisions? 
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The hydrodynamic (hydro) model applied to heavy ion data from the relativistic heavy ion col- 
lider (RHIC) in the form of single-particle spectra and correlations seems to indicate that a dense 
QCD medium nearly opaque to partons, a strongly-coupled quark-gluon plasma (sQGP), is formed 
in more-central Au-Au collisions, and that the sQGP may have a very small viscosity ("perfect liq- 
uid"). Measurements of radial and elliptic flows, with possible coalescence of "constituent quarks" 
to form hadrons, seem to support the conclusion. However, other measurements provide contradic- 
tory evidence. Unbiased angular correlations indicate that a large number of back-to-back jets from 
initial-state scattered partons with energies as low as 3 GeV survive as "minijet" hadron correla- 
tions even in central Au-Au collisions, suggesting near transparency. Two-component analysis of 
single-particle hadron spectra reveals a corresponding spectrum hard component (parton fragment 
distribution described by pQCD) which can masquerade as "radial flow" in some spectrum analysis. 
Reinterpretation of "elliptic flow" as a QCD scattering process resulting in fragmentation is also pos- 
sible. In this paper I review analysis methods and results in the context of two paradigms: the con- 
ventional hydrodynamics/hard-probes paradigm and an alternative quadrupole/minijets paradigm. 
Based on re-interpretation of fiducial data I argue that hydrodynamics may not be relevant to RHIC 
collisions. Collision evolution may be dominated by parton scattering and fragmentation, albeit the 
fragmentation process is strongly modified in more-central A-A collisions. 

PACS numbers: 12.38.Qk, 13.87.Fh, 25.75.Aq, 25.75.Bh, 25.75.-q, 25.75.Ld, 25.75.Nq, 25.75.Gz 



I. INTRODUCTION 

The hydrodynamic (hydro) model has been applied ex- 
tensively to heavy ion collisions at the super proton syn- 
chrotron (SPS) and RHIC as part of a search for for- 
mation of a quark-gluon plasma (QGP) [H, 0. Hydro is 
intended to describe A-A collision evolution in terms of 
flowing hot and dense matter, possibly a QGP Va- 
lidity of the hydro description could support inference 
of parton thermalization [J, |5| and direct comparisons 
with lattice QCD Q. The hydro model appears to be 
successful in representing some aspects of particle data. 
Hadron spectra have been described in terms of radial 
flow combined with the statistical model @, [!]• Some 
azimuth correlations have been described in terms of el- 
liptic flow [§]. Based on elliptic-flow systematics more- 
central RHIC Au-Au collisions have been characterized 
in terms of a strongly-coupled QGP (sQGP) with small 
viscosity — a "perfect liquid" [l(| EJ • 

According to the hydro model almost all particle pro- 
duction over much of A-A momentum space reflects rapid 
thermalization and development of a flow field in re- 
sponse to initial pressure gradients. Mass dependence of 
the hadron momentum distribution should reflect the un- 
derlying flow system [l2j. However, differential spectrum 
and correlation analysis reveals structures whose varia- 
tion with A-A centrality and collision energy contradict 
hydro expectations [H, EH EE EE E3) EH • Interpretation 
of measured V2 as the property of a flowing bulk medium 
can be strongly questioned 

In this paper I examine assumptions and procedures 
which seem to support a hydro description at RHIC and 
assess their validity in the context of fiducial data. I 
contrast the hydro description with a two-component 



model of elementary N-N collisions and the Glauber lin- 
ear superposition (GLS) reference model of A-A colli- 
sions. The emerging two-component phenomenology is 
consistent with QCD systematics, does not require the 
hydro model and conflicts with hydro in several ways. 

The paper is organized as follows: After sum- 
marizing two RHIC paradigms (hydro/hard-probes vs 
quadrupole/minijets) I review arguments for and against 
the hydro model. I summarize analysis methods and 
plotting formats which play a central role in shaping in- 
terpretation of RHIC data. I review the two-component 
model of nuclear collisions and the systematics of minijets 
(minimum-bias jets, Ej et ~ 3 GeV). I consider so-called 
"triggered" jet analysis and related systematic errors. I 
review assumed hydro initial conditions compared to re- 
sults from elementary and nuclear collisions. Lastly, I re- 
view hydro applications to single-particle spectra (blast- 
wave model vs fragmentation) and azimuth correlations 
(elliptic flow vs jet systematics). 



II. TWO RHIC PARADIGMS 

Two paradigms compete to describe RHIC data. The 
conventional RHIC paradigm emerged from the Bevalac 
heavy ion program, where nuclear collisions were domi- 
nated by semiclassical molecular dynamics, some degree 
of thermalization was achieved and thermodynamic state 
variables described collisions. Hydrodynamics plays a 
dominant role in the paradigm. Anticipated novel as- 
pects of RHIC collisions at much higher energies include 
possible formation of a thermalized partonic medium or 
QGP at larger energy densities and modification of par- 
ton hard scattering and jet formation by the medium. 
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An alternative paradigm is based on fragmentation 
processes in elementary e + -e~ (e-e) and N-N collisions as 
a QCD reference system. N-N collisions above y/s ~ 15 
GeV are described near mid-rapidity by a two-component 
model of longitudinal and transverse fragmentation. Lin- 
ear superposition of N-N binary collisions provides a ref- 
erence model for A-A collisions. Deviations from the ref- 
erence may reflect novel QCD physics in A-A collisions. 



A. The hydrodynamics/hard-probes paradigm 

In the conventional picture of RHIC collisions a ther- 
malized "bulk partonic medium" results from copious 
rescattering of ~1 GeV partons (gluons) 0, |2(J. Result- 
ing pressure gradients drive hydrodynamic expansion and 
formation of a velocity held (elliptic, radial and longitudi- 
nal flows). The medium expands, cools and "freezes out" 
to form hadrons which may continue to rescatter and ex- 
pand collectively. More-energetic partons (hard probes) 
are multiply scattered in the partonic medium and lose 
energy by gluon bremsstrahlung, possibly stopping in 
the medium (parton thermalization). The hadronic final 
state is described by hydrodynamic models, statistical- 
model hadrochemistry and perturbative QCD (pQCD) 
jet quenching [3|. 

Phenomenologically, hadron single-particle (SP) trans- 
verse momentum p t spectra (cf. Fig. 1 1 0[) are divided into 
a low-p t region < 2 GeV/c, nominally representing the 
thcrmalizcd bulk medium and described by blast-wave 
(BW) and statistical models, and a high-p t region > 5 
GeV/c dominated by parton scattering and fragmenta- 
tion and described by pQCD. The intermediate-p f re- 
gion (~2-5 GeV/c) manifests novel structure possibly 
described by hybrid models (e.g., recombination or co- 
alescence of "constituent quarks" ) denoted by ReCo [2l|, 

MM- 



B. The quadrupole/minijets paradigm 

In recent years differential analysis of spectrum and 
correlation data from p-p collisions has provided a de- 
tailed phenomenological reference for nuclear collisions. 
Accurate determination of A-A centrality extending to 
N-N collisions was also achieved. The combination es- 
tablished a Glauber N-N linear superposition or GLS ref- 
erence: What if nothing new happens in A-A collisions? 

Spectrum and correlation analysis of p-p collisions re- 
vealed a significant contribution from "minijets" (frag- 
ments from the minimum-bias parton spectrum domi- 
nated by ~3 GeV partons, cf. Sec. IVl]) contributing 
a significant fraction of the p-p hadron spectrum down 
to 0.3 GeV/c. Minijet phenomenology is fully consis- 
tent with jet systematics at larger energy scales — e + -e~ 
and p-p fragmentation functions [l4j. Minijet correla- 
tions demonstrate that observable consequences of QCD 



quanta at small energy scales persist in A-A collision data 
and even dominate the hadronic final state. 

The 2D angular autocorrelation method developed to 
study minijet correlations also determines (directly) the 
azimuth quadr upole moment (model-independent termi- 
nology (19 . Izil . [25|) conventionally identified with "el- 
liptic flow." The autocorrelation technique accurately 
resolves the azimuth quadrupole ("flow") from minijet 
structure ( "nonflow" ) . Quadrupole energy and centrality 
systematics appear to be inconsistent with hydro expec- 
tations or an equation of state (EoS) , instead are sugges- 
tive of QCD multipole radiation [251 ] . 

In the context of the GLS reference minijet and 
quadrupole variations with A-A centrality and energy 
challenge assumptions supporting the hydro model. 
Minijets serve as Brownian probes of hypothetical "bulk 
matter" [24| and test conclusions from the hydro model: 
Is there a collective flowing medium? Is the medium 
opaque to jets? Is viscosity relevant? This paper com- 
pares arguments supporting the two paradigms and chal- 
lenges the validity of a hydro description of RHIC colli- 
sions. 



III. ARGUMENTS RELATING TO HYDRO 

Application of hydrodynamics to RHIC collisions is in- 
herited from the Bevalac program which characterized 
nuclear collisions in terms of the molecular dynamics of 
nucleon clusters and nucleons and demonstrated collec- 
tive flow phenomena. Similar results were obtained at 
the AGS for systems of nucleons and nucleon resonances 
at the higher energy. Hydro arguments extended to the 
context of RHIC collisions are summarized below. 



A. Basic hydro arguments 

a. Initial conditions The pQCD energy spectrum 
for scattered partons is assigned a low cutoff energy 
(~ 1 GeV) — justified by saturation-scale arguments — 
which implies a large parton (gluon) phase-space density. 

b. Parton thermalization The resulting large parton 
density (10-30 times larger than expected from Glauber 
extrapolation of p-p collisions) implies thermalization 
through parton multiple scattering. 

c. Equation of state Assumed thermalization (=>• re- 
versibility, detailed balance?) justihes system description 
with state variables related by an EoS. 

d. Hydrodynamic evolution The EoS relates a large 
initial energy density to pressure gradients converted via 
hydrodynamic evolution to a bulk-matter velocity field 

e. Parton energy loss "Hard probes" of the ther- 
malized medium exhibit parton energy loss and modified 
fragmentation. A "well-understood" pQCD phenomenon 
is interpreted to reveal bulk-medium properties 
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/. Late hadronization The thermalized, flowing par- 
tonic bulk medium hadronizes by constituent-quark coa- 
lescence, producing nearly all final-state hadrons. 

g. Hadron rescattering Isentropic transverse expan- 
sion of the resulting hadron resonance gas leads to differ- 
ent chemical and kinetic decoupling temperatures. 

h. Final-state flows Evidence is sought in final-state 
hadronic spectra and correlations — e.g., mass depen- 
dence of (pt), radial and elliptic flows — for a flowing par- 
tonic medium established prior to hadronization. 



B. Contrasting arguments 

a. Parton spectrum A parton spectrum cutoff near 
3 GeV is established by minijet correlations and spectrum 
hard components from p-p and Au-Au collisions. 

b. Thermalization Three separate components [soft, 
hard (minijets) and quadrupole] have distinct character- 
istics, different mechanisms and remain largely indepen- 
dent throughout the collision for all A-A centralities. 

c. State variables There is no evidence for a ther- 
modynamic state (homogeneity, reversibility or detailed 
balance) and therefore no support for state variables. 

d. Parton multiple scattering Minijet correlations 
indicate that essentially all initial-state scattered par- 
tons appear as correlated structures in the hadronic final 
state, even in central Au-Au collisions. There is negli- 
gible parton rescattering — minijet azimuth widths grow 
smaller, not larger with increasing A-A centrality. 

e. Flows No radial flow is observed in p t spectra. 
SP spectrum structure is described by two fragmentation 
components, is inconsistent with the blast-wave model. 
i>2 (Pt) data can be interpreted in terms of QCD multipole 
radiation. 

In the context of the two-component model, arguments 
which support hydro implicitly invoke a competition be- 
tween soft (nucleon fragmentation) and hard (scattered- 
parton fragmentation) components for final-state hadron 
production near mid-rapidity. The spectrum soft com- 
ponent, the pt spectrum from projectile nucleon frag- 
mentation, represents the majority of hadron production 
for all A-A centralities. It is a universal feature of frag- 
mentation observed also in e-e collisions. The spectrum 
hard component, the pt or y t distribution of fragments 
from minimum-bias large-angle-scattered partons, repre- 
sents the main transport mechanism from longitudinal to 
transverse parton phase space. 

Arguments supporting hydro assume that most of 
the hard component is thermalized, and the resulting 
medium dominates hadron production in more-central 
collisions. Collisions are in some sense "opaque" to scat- 
tered partons which undergo many rescatterings. The 
N-N soft component likewise dissolves into the medium 
in more-central Au-Au collisions. Hadron production is 
dominated by hadronization of the thermalized medium. 
A small jet-correlated fraction results from fragmenta- 
tion of rare high-p t partons (possibly a "surface bias" ) . 



Observed strong minijet correlations and independent 
quadrupole systematics contradict those assumptions. 

IV. ANALYSIS METHODS 

Very different interpretations can emerge from the 
same data depending on analysis methods and plot- 
ting formats. Preferred analysis methods extract all 
information from spectra and correlations in directly- 
interpretable form. Differential analysis incorporates 
well-understood and comprehensive reference models of 
nucleon and parton scattering and fragmentation, A-A 
centrality dependence and two-particle correlations. 

A. References 

a. Glauber model A-A collision geometry parame- 
ters n part , Ubin, b and cr/cro related to observables such 
as n c h at mid-rapidity; the preferred centrality measure 
is mean participant pathlength v = 2ni, in /n part (2r| . 

b. Parton scattering and fragmentation Complete 
parametrization of parton scattering and fragmentation 
for e-e and p-p collisions; parton power-law spectrum in- 
ferred from p-p data; calculated fragment distributions 
derived from a pQCD folding integral fbH. l27l]. 

c. Two-component model Soft+hard spectrum and 
correlation components representing non-single diffrac- 
tive (NSD) N-N soft scattering and (longitudinal) frag- 
mentation and pQCD semihard parton scattering and 
(transverse) fragmentation [l3], Ha HH ■ 

d. Glauber linear superposition (GLS) Essential ref- 
erence for centrality dependence of A-A collisions over 
the complete centrality range; spectrum and correla- 
tion measures normalized by participant-pair number 
npart/2, reference soft component is invariant, refer- 
ence hard component increases proportional to nunary'i 
physics unique to A-A deviates from the GLS reference. 

B. Plotting formats 

In high-energy collisions the choice of plotting 
variables — yt vspt, y vs x v or £ p , v vs n part , V2 vs Ap[2]— 
strongly affects the visibility of important structure and 
its physical interpretation. Conventional momentum- 
space variables are y z for longitudinal momentum and 
p t for transverse momentum. A typical mid-rapidity de- 
tector acceptance at RHIC (STAR TPC) corresponds 
(for pions) to \p z \ < 0.2 GeV/c (\y z \ < 1) and p t < 10 
GeV/c (yt < 5). A more sensible choice might be y z 
and yt or even p z and yt. Linear pt over-emphasizes 
the "high-p t " region where physical changes are modest 
at the expense of the sm&\\-p t region where substantial 
new QCD physics emerges at RHIC. Transverse rapidity 
y t = log{(p t + m t )/m } ~ \og(p t ) more-clearly displays 
systematic QCD trends. 
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Two examples of variable choices and consequences are 
presented, one from parton fragmentation, the other from 
differential spectrum analysis with ratios. 

a. Parton fragmentation The most important kine- 
matic region for parton fragmentation in nuclear colli- 
sions includes small parton energies and small fragment 
momenta where most fragments are produced. Fragmen- 
tation functions (FFs) are conventionally represented on 
linear momentum fraction x p — p/pjet or £p = ln(l/x p ), 
whereas the format most relevant to nuclear collisions is 
on rapidity y — ln{2(£' + p)/mo} or normalized rapidity 
u = (y-y m in)/(ymax-y m in), with y m „ x = ln(2p jet /m ) 
and y m i n ~ 1/3 for e + -e~ collisions [27[. 
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FIG. 1: Left panel: Beta distribution (solid) [27| and KKP FF 
(dashed) [3(j curves compared to OPAL 91 GeV data points 
(open circles) [2|| on linear momentum variable x p . Right 
panel: The same curves and data transformed to normalized 
rapidity u. The vertical dotted lines correspond to x v — 0.1. 

Fig. Q] illustrates the consequences. Fiducial FF data 
from OPAL (y/s = 91 GeV) [H[ are plotted on linear 
momentum fraction x p on the left, and on normalized 
rapidity u on the right. Less than 10% of the distribu- 
tion falls to the right of the dotted line in each panel. The 
dashed curves represent a 15-parameter pQCD-theory 
formulation on x p [30(. The solid curves describe a two- 
parameter beta distribution on u [27| . While more-recent 
theory describes FF data down to x p ~ 0.05 [3l| this fig- 
ure illustrates visual consequences of plotting formats. 

In the left panel pQCD (KKP) and beta parameteriza- 
tions appear to describe the data equally well for x p > 0.1 
where pQCD DGLAP evolution is emphasized [271] . How- 
ever, the relation between models and data for 90% of the 
fragments, to the left of the dotted line (the region which 
dominates nuclear collisions), is not resolved visually. In 
the right panel the large discrepancy between theory and 
data becomes clear, as does the accuracy of the beta 
parametrization which describes all fragments over all en- 
ergy scales relevant to nuclear collisions (p t > 0.1 GeV/c, 
Ejet > 3 GeV) [13]. FF modes occur in the interval 1-2 
GeV/c for all parton energies up to several hundred GeV. 
In nuclear collisions at least 50% of scattered-parton frag- 
ments fall below 1 GeV/c. 

b. Spectrum ratios Spectrum ratios are differential 
measures of spectrum centrality variation. Two ratio 



measures have been defined 
1 Paa 



Raa — 
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where v = 2nbi n /n part is the participant-nucleon mean 
pathlength, Snn and Hnn are reference per-participant- 
pair soft and hard spectrum components respectively, 
and Haa is the hard component extracted from data [13j . 

Fig. [2] (left panel) shows ratio Raa for five centrali- 
zes of Au-Au collisions (five bold curves of different line 
styles) plotted on transverse rapidity y t [13j . Also shown 
are five reference curves (light solid curves) with Haa 
replaced by reference Hnn in Eq. (p} (first line), and the 
solid points with Haa replaced by H pp data from p-p 
collisions [35| corresponding to v = 1 which fall along 
the reference line at 1. 
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FIG. 2: Left panel: Spectrum ratio Raa for five centralities 
of 200 GeV Au-Au collisions (bold curves) and 200 GeV p-p 
collisions (points) plotted on transverse rapidity yt The 
thin solid curves are Glauber linear-superposition (GLS) ref- 
erences, v is a centrality measure. Right panel: Alternative 
spectrum ratio taa for the same data. The GLS reference 
for all cases is unity. The structure at smaller yt is newly 
revealed. 



Fig. [3 (right panel) shows alternative ratio taa de- 
fined in terms of spectrum hard components from the 
same data as in the left panel (l3j ]. Soft component 
Snn is eliminated from the ratio, and full access to hard- 
component centrality trends is thus achieved. The cor- 
rect reference is unity for all yt, making deviations from 
the reference unambiguous. More importantly, sensitiv- 
ity to deviations is uniform over all yt- The 200 GeV 
p-p and 60-80% central Au-Au data clearly agree with 
the reference. The large excursions at smaller y t for 
more-central Au-Au collisions, apparent for the first time 
with rAA-, represent most of the 30% increase in per- 
participant multiplicity for more-central collisions. 

Fig. [5] and Eq. demonstrate that Raa can be 
quite misleading. The soft component included in Raa 
strongly suppresses hard-component structure below 4 
GeV/c, and the true reference curves are not acknowl- 
edged (the reference is usually assumed to be unity). At 
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Vt ~ 2 (pion pt ~ 0.5 GeV/c) 

1 , Hnn 
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where H^n / Snn ~ 1/170 for pions, and con- 
tains all information on spectrum centrality evolution. 
Hard-component centrality dependence below 2 GeV/c 
(yt ~ 3.3) is therefore strongly suppressed visually, and 
the relation to the correct reference is obscured. Plot- 
ting Raa on conventional pt further distracts from rele- 
vant structure at smaller y t , emphasizing high-p t "hard 
probes" to the exclusion of important new fragmentation 
structure in a pt interval nominally assigned to hydro. 

The upper-most (solid) data curve (60-80% central) 
in the left panel could be described as "suppressed," 
but is actually consistent with the GLS reference (up- 
permost thin solid curve) over the entire yt acceptance. 
The taa trend (right panel) for more-central collisions 
is still "suppression" above yt = 4 (pt = 4 GeV/c), but 
important new large-amplitude structure appears at the 
left end of the spectra. Centrality dependence near 10 
GeV/c (yt = 5, suppression) is closely related to central- 
ity dependence below 0.5 GeV/c (yt — 2, enhancement). 
No spectrum structure in this differential analysis corre- 
sponds to radial flow [l3| (and cf. Sec. IIX Cj) . 



C. Hydro-motivated analysis methods 

Figure [3] illustrates variable choices and plotting for- 
mats critical to hydro interpretations. The same V2(pt) 
data for three hadron species from minimum-bias Au-Au 
collisions at 200 GeV appear in each panel [H, [33| . In 
the left panel is the conventional «2(pt) vs pt format. In 
the right panel the data have been processed so as to 
reveal quadrupole spectra on transverse rapidity y t for 
each hadron type. 

V2(pt) is proportional to the ratio of two spectra: 
quadrupole spectrum P2(yt] Ay t $, T%) in the numerator 
and the single-particle or SP spectrum in the denomi- 
nator. The quadrupole spectrum can be inferred from 
measured quantities by [19( 



p(yt 



V2{Vt 
Pt 



Pt 



Ptlt(l ~ Pt) 

f(yf, Ay i0) Ay t2 



7t(l~A) 
Ay t 2 p 2 (yu Ay m ,T 2 



(3) 



p(Vt) = Po(yt) + P2(Vu <t>) is the SP spectrum appearing 
in the denominator of V2(yt)- Monopole boost Ay to (in- 
ferred from data) is the transverse boost of a common 
source for all hadrons associated with the quadrupole. 
The quantities in curly brackets are determined by Ayto 
[except T2 which is inferred from the shape of P2(Vt)\- 
Function f(yt) is unity for smaller y t but increases 
smoothly to about 1.2 at yt ~ 3 for the plotted data [lj|. 

Fig. [3] (left panel) shows the conventional plotting for- 
mat for V2(pt)- The mass ordering at smaller p t is inter- 
preted to imply hydro expansion. "Saturation" at larger 
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FIG. 3: Left panel: pt-differential elliptic flow measure V2(pt) 
for three hadron species from minimum-bias 200 GeV Au- 
Au plotted in a conventional format [33, Hydro the- 
ory curves A and B are described in the text. Right panel: 
The same «2(pt) data plotted in an alternative format using 
proper rapidity yt for each mass species, which reveals az- 
imuth quadrupole spectra with Levy form and common trans- 
verse boost Ayto ~ 0.6. 



Pt might indicate jets correlated with the reaction plane. 
Sec. IX Bl presents a more-detailed discussion of V2(pt)- 
The full SP pt spectrum in the i>2 denominator contains 
contributions from parton scattering and fragmentation 
(hard component) which dominate V2(pt) at larger p t . 
Theory curves A and B are described in Sec. IX Bl 

Fig. 02 (right panel) based on Eq. removes from 
v 2 (pt) the SP spectrum in the denominator as well as 
an extraneous kinematic factor p t to reveal quadrupole 
spectra. Plotting data on rapidity proper to each hadron 
species reveals a common boosted source (aligned left 
edges) and source boost Ay t Q. The universal Levy form 
of the spectra (solid curves) is also apparent, permitting 
direct interpretation of V2(pt) data without the interme- 
diary of the theory to be tested. 



D. Supporting material 



Support for the arguments in this paper falls into sev- 
eral categories. Basic analysis methods include central- 
ity determination (26[, 2D angular autocorrelations and 
inversion of fluctuation scale dependence 34] . Determi- 
nation of the two-component structure of single-particle 
(SP) spectra is described in [H, [35|. Fragmentation in 
elementary and A- A collisions is discussed in [3, l27j |. 
Methods for two-particle correlations in p-p collisions are 
presented in (28l. |36|. Correlation methods for A- A col- 
lisions are presented in [ll| [lj], HH • Analysis of trig- 
gered azimuth correlations is discussed in [33. Analysis 
methods related to the azimuth quadrupole (elliptic flow) 
are presented in 0, [H, HI . 
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V. THE TWO-COMPONENT MODEL 

The two-component (soft+hard) model of hadron pro- 
duction is an essential reference for spectra and corre- 
lations (setting aside the azimuth quadrupole — "elliptic 
flow" — as an independent and relatively small third 
component). The components represent 1) non-single 
diffractive (NSD) N-N soft scattering and (longitudi- 
nal) fragmentation and 2) large-angle pQCD semihard 
parton scattering and (transverse) fragmentation respec- 
tively. Coupled to Glauber linear superposition the two- 
component model provides a comprehensive reference for 
A- A collisions. 



A. Spectrum model for p-p and A-A collisions 

The two-component spectrum model for p-p collisions 
sorted according to detected multiplicity h c h ~ n c h/2 is 



1 1 dn ch 



n s (n ch ) y t dy t 



S (y t ) + ^p^H (y t ), (4) 
n s {n ch ) 



with rih/n s rs ah c h and a ~ 0.01. Unit-normal model 
functions So and Ho are extracted from data in [3||. 
Spectrum components Snn — n s Sq and H n n = rihHo 
((^-integrated 2D densities) describe per-participant-pair 
particle densities (n x ~ dn x /drj are defined in one unit 
of rf). Soft component So(yt) is derived from spectrum 
data as the physical-model-independent limit 



c , < 1 1 dn ch 

&o{yt) = hm — . 

n c h Vt dy t 

Hard component Hq is defined by the average of 



(5) 



a rich H (y t ) 



1 



1 dn ch (h ch ) 



n s (n ch ) y t dy t 



Sofa), (6) 



with n s = n ch /(l + ah ch ). 

Figure [4] (left panel) shows pt spectra for ten multiplic- 
ity classes from NSD p-p collisions at y/s = 200 GeV [35| • 
The spectra are normalized to soft-component multiplic- 
ity n s determined by an iterative process. So is clearly 
the limit of measured spectra as h c h — > 0, per Eq. ([5]). 

Fig. 2] (right panel) shows the two-component model 
of Eq. 6D with a = 0.01 and S a and H from |f|. The 
n c /i-independent model functions represent all significant 
structure in the spectrum data. The two-component 
model is a nontrivial decomposition of p-p pt spectra. 

The Glauber model of A-A collisions relates mean 
participant-pair number n part /2 and mean number of bi- 
nary N-N collisions nu n to the fraction of total cross 
section ct/cto or impact parameter b (26[. The fractional 
cross section in turn relates the Glauber parameters to 
an observable such as multiplicity n c h- Mean partici- 
pant pathlength v = 2nbi n /n pa rt is the prefered central- 
ity measure for the GLS reference. The two-component 




FIG. 4: Left panel: Single-particle (SP) spectra from 200 
GeV p-p collisions for ten multiplicity classes normalized by 
soft-component multiplicity n s j35|]. Unit-normal Levy distri- 
bution So is the soft-component model function common to 
all spectra. Right panel: Two-component model for spectrum 
data in the left panel. Unit-normal hard component Ho is a 
Gaussian, with exponential tail representing a pQCD power 
law on pt. 



model for A-A spectra is then formulated by analogy with 
p-p h c h dependence 



1 dn r 



n part yt dy t 



SNN{yt) + vHAA{yt,v), (7) 



where the soft component is by hypothesis indepen- 
dent of centrality, and all centrality dependence devi- 
ating from the GLS is absorbed into hard component 
Haa{Vu v ) [uj - The Glauber linear-superposition (GLS) 
model corresponds to H aa — > Hnn in Eq. ([7]). 



B. Spectrum soft component 

Fig. E] (left panel) shows a pion spectrum (2 x ir~) on 
nit (solid points) for SPS fixed-target 0-2% central S-S 
collisions at 200A GeV (y/s^N = 19.4 GeV) [38]. A hy- 
dro model applied to those data was used to infer radial 
flow with mean transverse speed (Pt) ~ 0.25 [7|. The 
curve labeled Snn is the soft component for yfs = 200 
GeV NSD p-p collisions (from Sq in Fig. [4j transformed 
to m t . The line labeled M-B is the Maxwell-Boltzmann 
limiting case of Snn- For comparison, a n + + n~ spec- 
trum from p-p collisions at 158A GeV (^/s — 17.3 GeV) 
is included (open circles) [3^| . Small differences at larger 
m t may be due to a combination of initial-state k t effects 
and semihard parton scattering [40]. Is there significant 
radial flow in p-p collisions, or is the hydro model not 
relevant to S-S data? 

Fig. [5] (right panel) shows an mt spectrum (solid dots) 
from LEP e+-e~ (e-e) collisions at y/s = 91 GeV [H|. 
The spectrum is derived from a sphericity analysis of q-q 
dijets and estimates the fragment momentum distribu- 
tion transverse to the q-q axis. LEP dn c h/dpt data have 
been normalized as indicated to match p-p soft compo- 
nent Snn (dashed curve), and hence the SPS S-S spec- 
trum. The shape of the wit spectrum from 91 GeV LEP 
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FIG. 5: Left panel: mt spectra from 0-2% central S-S colli- 
sions at 200A GeV (solid points, -^/sjvjv = 19.4 GeVn3a] and 
from p-p collisions at ^/s = 17.3 GeV (open circles) [3f|. The 
dashed curve is the Levy soft component from 200 GeV p-p 
collisions. The dash-dotted curve is the Maxwell-Boltzmann 
limiting case of 5jvjv- Right panel: mt spectrum (points) from 
e + -e~ collisions at yfs = 91 GeV piU ]. Curves are duplicated 
from the left panel. The data are normalized as indicated in 
the axis label to match the Snn curve at larger mt. 



FIG. 6: Left panel: Hard components (solid curves and 
points) from ten multiplicity classes of 200 GeV p-p colli- 
sions [35|]. The data are normalized to non-single diffractive 
(NSD) p-p collisions (observed h c h ~ 1.25 in one unit of rf). 
The dash-dotted curve is Gaussian+exponential tail reference 
Ho- Right panel: p-p hard component (points) averaged over 
several h c h classes and parton fragment distribution (FD, 
solid curve) obtained from a pQCD folding integral. Dot- 
ted curves correspond to ±10% shifts in the ~ 3 GeV parton 
spectrum endpoint. 



FFs is thus consistent with the soft component of p-p 
spectra at 200 GeV and with the full spectrum from cen- 
tral S-S collisions at 19 GeV. 

The commonality of the soft-component spectrum 
shape (Levy distribution) across energies and collision 
systems suggests that the soft component is a universal 
feature of fragmentation whatever the leading particle — 
parton or hadron. It is unlikely that hydro expansion 
plays a role in LEP e + -e~ collisions. Thus, inference of a 
radial flow velocity from the same (Levy) spectrum shape 
in A-A collisions is questionable. 



C. Spectrum hard component 

Figure [5] (left panel) shows spectrum hard-component 
data from 200 GeV NSD p-p collisions in the 
form [nh(1.25)/nh(n c h)] H pp for ten multiplicity h c h 
classes [35[. The normalization corresponds to NSD 
value n c h — 2.5 in one unit of r\. The common 
shape, well-described by a Gaussian with exponential tail 
(dash-dotted curve), is plotted as H in Fig. [4] (right 
panel). The exponential on y t represents QCD power 

law Pt - n<3CD [rim. 

Figure [6] (right panel) shows p-p hard-component data 
averaged over several multiplicity classes (solid points) 
compared to a pQCD fragment distribution (FD, solid 
curve) and the Gaussian-plus-tail reference (dash-dotted 
curve). The FD was calculated by combining p-p frag- 
mentation functions [42j with a power-law parton spec- 
trum (Fig. O left panel, solid curve) in pQCD folding in- 
tegral Eq. [14| . This comparison confirms that the p-p 
spectrum hard component is a minimum-bias fragment 
distribution (minijets). The hadron spectrum data de- 
termine the parton spectrum lower-cutoff energy to 5%. 
In A-A collisions hard component Haa evolves dramati- 



cally, revealing suppression at larg er p t and much larger 
enhancement at smaller p t [H, H3| (and cf. Fig. [T0|) . 



D. Two-component correlations 

The soft+hard two-component decomposition applies 
also to two-particle correlations on angular difference 
variables (t]a, 4>a) ( e -S-; ?7A = Jli — V?) an( i transverse 



rapidity (y t ,Vt) 



In p-p collisions the 



soft component is confined to unlike-sign pairs peaked 
at the origin on r)& and nearly back-to-back on azimuth, 
the structure being consistent with longitudinal (projec- 
tile nucleon) fragmentation to charged-hadron pairs. 

The hard-component correlation structure has the 
properties expected for (mini)jets: a same-side (SS, 
I'M | < 7r/2) 2D peak at the angular origin and an away- 
side (AS, |0a| > t/2) ridge corresponding to back-to- 
back (0a = 7r) parton scattering. On (yt,yt) the hard- 
component peak mode at (2.7,2.7) (p t ~ 1 GeV/c) is con- 
sistent with the single-particle spectrum hard component 
in Fig. [6l p-p structure agrees semi-quantitatively with 
PYTHIA correlations [Hli|. 

In Au-Au collisions the correlation soft component 
decreases rapidly to zero with increasing collision cen- 
trality. Hard-component correlations change dramati- 
cally: Both angular and (yt,Vt) correlations follow the 
GLS reference up to a particular centrality, then de- 
velop large deviations from GLS (increases) for more- 
central Au-Au collisions [n| which correspond to a 
sharp transition at the same centrality for the spec- 
trum hard component [H, Hi]. The increased jet struc- 
ture contradicts triggered dihadron correlations which 
imply strong jet s upp ression (parton thermalization) (cf. 

Sec. En} jli, EIIII 
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VI. MINUETS: MINIMUM-BIAS JETS 

The term "minijet" refers experimentally to structure 
in hadron spectra and correlations produced by frag- 
mentation from the minimum-bias energy spectrum of 
small- a; partons (mainly gluons) scattering to large an- 
gles near mid-rapidity. No pt condition is imposed on de- 
tected hadron fragments. The inferred parton spectrum 
is dominated by partons near 3 GeV (Q ~ 6 GeV), each 
parton on average fragmenting to two charged hadrons 
with p t ~1 GeV/c. Parton spectrum structure and corre- 
sponding fragment distributions (FDs) in hadron spectra 
and correlations are described in 1141. 



A. Minijets and theory 

The minijet concept originated with an analysis of the 
spectrum of event-wise Et clusters down to 5 GeV in 
calorimeter data from 200 GeV p-p collisions [13] • pQCD 
analysis of the UA1 (mini) jet spectrum assumed a par- 
ton spectrum cutoff near 3 GeV to obtain a jet total cross 
section of 2-3 mb. flR |4!^]. A lower estimate of the cut- 
off for RHIC A- A collisions (1 GeV) has been based on 
saturation-scale-model (SSM) arguments [l.[20|. 

Minijets were assumed by some theorists to be un- 
resolvable in the A-A final state due to thermaliza- 
tion, but should contribute a substantial fraction of 
hadron and pt /E t production in more-central Au-Au col- 
lisions 0, Ho, H] . In a hydro theory context thermalizcd 
minijets are expected to provide the initial energy density 
required to drive hydrodynamic expansion In effect, 
almost all particle, pt and E t production in the final state 
of central Au-Au collisions is attributed to thermalized 
minijets. 

Other descriptions of minijets were intended to distin- 
guish true QGP manifestations in RHIC collisions from 
pQCD processes. Minijets were modeled with the HI- 
JING Monte Carlo [13, [Hj. The default HIJING parton 
spectrum cutoff is po = 2 GeV/c. HIJING without "jet 
quenching" is a Glauber linear-superposition (GLS) ref- 
erence which provides a semi-quantitative description of 
observed minijet phenomenology in p-p and peripheral 
Au-Au collisions down to pt ~ 0.1 GeV/c hadron mo- 
mentum. 



B. Minijets and RHIC data 

Systematic studies of minijets in p-p and Au-Au colli- 
sions at RHIC [TH, [l6|, [H, Hy] are inconsistent with the- 
oretical assumptions of minijet thermalization. Spectra 
and correlations provide substantial evidence that back- 
to-back minijets are actually resolved in the final state of 
A-A collisions, are not "thermalized" even in central Au- 
Au collisions, although parton scattering and fragmen- 
tation are strongly modified. Novel techniques leading 



to that conclusion include accurate centrality measure- 
ment from N-N to b — Au-Au [26| . differential two- 
component analysis of p-p and A-A spectra 131 35[ and 
2D angular autocorrelations [H, [13, El S [HJ3J 

Application of the two-component spectrum model 
to Au-Au collisions reveals soft+hard structure simi- 
lar to p-p collisions [l3|. The hard-component abun- 
dance relative to participant pairs increases at least as 
fast as oc v = 2 n binary / n par u cipant (mean participant 
path length) representing the N-N binary-collision scal- 
ing expected for parton scattering and fragmentation. 
Hard-component centrality dependence near 0.5 GeV/c 
(yt ~ 2) (large increase) closely corresponds to that near 
10 GeV/c (y t ~ 5) (smaller decrease) associated with 
parton energy loss or "jet quenching," consistent with the 
hard component being a minimum-bias parton fragment 
distribution extending down as far as 0.1 GeV/c fl4j ]. 

Angular correlations identified as minijets increase dra- 
matically in more-central Au-Au collisions, in contrast 
to strong "jet quenching" inferred from Raa and par- 
ton thermalization inferred from triggered dihadron cor- 
relations. The close correspondence of trends for minijet 
correlations and spectrum hard component is reasonable 
since SP spectra include the marginal projections of jet 
correlations. 

Observed minijets reveal the underlying parton spec- 
trum and the "initial state" of A-A collisions. Compar- 
ison of spectrum hard components with calculated FDs 
indicates that the parton spectrum cutoff is much higher 
(3 GeV) than what is inferred from saturation-scale ar- 
guments that accommodate hydro (1 GeV). Survival of 
almost all minijets to hadronic decoupling contradicts hy- 
dro expectations. 

The GLS reference and ideal hydro in an opaque 
medium are limiting cases of A-A collisions which de- 
fine a metric between perfect transparency and complete 
opacity to partons. Where do A-A collisions of given 
centrality and energy lie on that axis? What is the quan- 
titative deviation from GLS transparency? Minijets act 
as "Brownian probes" to answer such questions [2J| . 



VII. "TRIGGERED" JET ANALYSIS 

As an alternative to event-wise jet reconstruc- 
tion in the high-multiplicity A-A environment trig- 
ger/associated pt cuts are imposed on azimuth correla- 
tions to visualize jets combinatorially. The cuts bias the 
underlying parton energy spectrum, the SP spectrum of 
hadron fragments and their angular correlations. 

Analysis and interpretation of triggered hadron corre- 
lations is challenging. Conventional trigger techni ques 
result in substantial underestimation of jet yields [37l. |44|. 
l45l |4(| and inference of a dense medium formed in more- 
central collisions [H, [HI] , with parton multiple scatter- 
ing, energy loss and production of Mach shocks in the 
medium [54. l55l l56l]. 

Figure [7] illustrates basic issues for triggered-jet data 
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analysis. Pairs of particles are combined from "trigger" 
and "associated" p t bins, the trigger p t bin being higher 
and usually disjoint from the associated bin. The number 
of pairs is normalized by the number of trigger particles. 
The basic pair distribution then contains a large fraction 
of uncorrelated combinatoric pairs and a small fraction 
of correlated pairs from two primary sources: azimuth 
quadrupole ("elliptic flow" measured by V2) and jets. 




FIG. 7: Left panel: Simulated "raw" (unsubtracted) dihadron 
correlations (points), azimuth quadrupole sinusoid with am- 
plitude A% (light dotted curve), and corresponding ZYAM 
subtracted background (bold dotted curve). The ZYAM- 
estimated background offset is the dotted line. Right panel: 
Result of ZYAM background subtraction in the left panel 
(points and bold curve). Dash-dotted and dashed curves are 
same-side (SS) and away-side (AS) jet peaks input to the 
simulation. The minimum at n in the AS (away-side) peak is 
notable. 

Conventional analysis combines independent v 2 mea- 
surements with the ZYAM (zero yield at minimum) con- 
vention to separate combinatoric background from jet 
correlations. In Fig. [JJ (left panel) the light dotted 
curve illustrates the azimuth quadrupole amplitude in- 
ferred by fitting the entire distribution with the form 
A + A 2 cos(2A</>), with A 2 /2A = uf{2} defining v 2 {2}. 
For more-central collisions where the non-jet quadrupole 
amplitude is small v|{2} 0.2 Ass /2, where Ass is 
the amplitude of the same-side (Gaussian) jet peak [37j |. 
Jet correlations then dominate v 2 {2}. The combination 
of «2{2} and ZYAM defines the bold dotted curve sub- 
tracted from the "raw" pair distribution to estimate jet 
structure. 

Fig. [7J (right panel) shows typical results (solid curve 
and points) for more-central A- A collisions. The dash- 
dotted and dashed curves show the SS and AS jet peaks 
input to the simulation. The dotted sinusoid is the ex- 
traneous quadrupole component (jet-induced v 2 {2}) im- 
posed in the ZYAM subtraction. The dotted line indi- 
cates the correct offset. From this simulation it is clear 
that the conventional background subtraction procedure 
can reduce the true jet yield by a large fraction. The AS 
peak, distorted by v 2 oversubtraction, acquires a mini- 
mum at 7T interpreted in terms of Mach shocks or cones. 

Fig . [5] (left panel) shows data from Fig. 1 (upper left) 
of 0. ZYAM-subtracted data from 200 GeV 0-12% cen- 
tral Au-Au collisions for (trigger x associated) 4-6x0.15-4 



GeV/c pt cuts are shown as solid points. Correspond- 
ing p-p data are shown as open circles. The bold solid 
curve is a free fit to the Au-Au data with offset (PI), 
SS Gaussian (amplitude P2, width P3), AS dipole (P4) 
and quadrupole (P5). The resulting offset is the solid 
line at Pi, the SS Gaussian is the dash-dotted curve, the 
dipole is the dashed curve, and the (negative) quadrupole 
is the dotted curve. The impression can be formed that 
jet yields in central Au-Au collisions are comparable to 
those in p-p collisions because of strong jet quenching in 
a dense medium. 
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FIG. 8: Left panel: ZYAM-subtracted angular correlations 
for 0-12% central 200 GeV Au-Au collisions (solid points) 0] 
with free fit (bold solid curve) of SS Gaussian (dash-dotted 
curve), AS dipole (dashed curve) and quadrupole (dotted si- 
nusoid). Open points are p-p data relative to ZYAM zero. 
Right panel: ZYAM subtraction reversed, true zero level re- 
covered from free fits to data (solid points), compared to p-p 
data (open symbols). 

Fig. |U (right panel) shows a reconstruction of the origi- 
nal ( "raw" ) pair distribution prior to ZYAM subtraction 
based on Pi and P5 from the left panel. The relation of 
p-p to central Au-Au data is quite different. Both SS and 
AS peaks increase by a factor six from p-p to central Au- 
Au, a large increase in the jet yield which is not apparent 
in the left panel determined by conventional ZYAM sub- 
traction. Strong suppression of SP hadron spectra at 
larger p t is apparently more than compensated by en- 
hancement of jet yields at smaller p t . 



VIII. HYDRO AND A-A INITIAL CONDITIONS 

Calculations of hydrodynamic evolution in nuclear col- 
lisions must specify the "initial conditions" — the thermo- 
dynamic state at some initial time, including the energy 
density, velocity and matter density distributions [3[ . For 
hydrodynamics to be a valid description of nuclear col- 
lisions a large number of initial-state scattered partons 
must thermalize to produce the energy density required 
to drive subsequent hydro evolution. Parton spectra ex- 
tending down to 1 GeV and copious rescattering are re- 
quired. pQCD can predict the shape of the initial parton 
energy spectrum, but the effective lower limit or cutoff of 
the spectrum, and the extent of subsequent parton ther- 
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malization, must be inferred from the hadronic final state 
and/or supplementary theoretical arguments. 



A. Initial conditions from theory 

Event-wise analysis of UAl Et distributions from 200 
GeV p-p collisions revealed minijet structure down to 5 
GeV jet energy, later interpreted from calorimeter back- 
ground estimates as due to partons down to 3 GeV [47l ]. 
pQCD analyses which assumed a minijet spectrum cut- 
off near 3 GeV estimated a jet cross section ~ 2 — 4 mb. 
[48l |49| . A 3-GeV cutoff is consistent with the observed 
p-p p t spectrum hard-component at 200 GeV [l4[ ■ 

Estimates of the spectrum cutoff for RHIC A-A colli- 
sions have been based on saturation-scale model (SSM) 
arguments 0, [20] | . Given a power-law form of the parton 
spectrum nearly all scattered partons emerge at the cut- 
off energy. Parton-related physics is then dominated by 
that (saturation-scale) energy. The SSM argument main- 
tains that parton scattering is limited only by the initial 
flux of partons from nucleon/nucleus projectiles. The 
parton spectrum cutoff is then determined by (gluon) 
saturation in the parton distribution functions (PDFs) of 
the colliding nuclei (or nucleons). SSM arguments lead 
to cutoff energy 1 GeV at ^/sWn = 200 GeV 0]. 

The SSM-based cutoff implies large parton phase-space 
densities which could be consistent with secondary par- 
ton scattering and partial thermalization (20l |. But the 
predicted hadron multiplicity and pt systematics are in- 
consistent with measured hadron spectra and correla- 
tions. And, the SSM argument that only initial-state 
saturation limits parton scattering may be incomplete. 
Hadronization (hadron density of final states) may be 
the determining factor in the parton spectrum cutoff, as 
observed in p-p collisions [l4T |. 

Figure [9] (left panel) shows parton power-law spectra 
with cutoff near 3 GeV (y m ax = l n (2iy e t/m 7r ) = 3.75) 
derived from comparisons with p-p and Au-Au spectra 
(solid and dash-dotted curves respectively) [l4j]. The 
bold dotted curve is an ab initio pQCD calculation ex- 
tending down to 1 GeV [20| . The absolute magnitudes 
agree well near 3 GeV, but large differences in parton 
and hadron yields arise from the different cutoffs. 



B. Initial conditions from experiment 

An alternative estimation of initial conditions can be 
obtained by direct comparison of calculated fragment dis- 
tributions (FDs) with spectrum hard components. The 
FD is defined by the pQCD folding integral [3] 



d 2 n h _ e(6rj,Arj) 
dy drj o- NS d At] 



dymax ^xx(|/i yn 



da, 



dijet 



dy„ 



,(8) 




FIG. 9: Left: Dijet (parton-pair) transverse energy spectra 
on rapidity ymax ~ hi(2 Ej e t/m n ). The solid curve is deter- 
mined by a measured p-p spectrum hard component. The 
dash-dotted curve illustrates reduction of the cutoff energy 
inferred for central Au-Au collisions. The bold-dotted theory 
curve labeled pQCD was derived from [2(j. The light dot- 
ted extrapolation to 1 GeV corresponds to a saturation-scale 
cutoff estimate 0, [2(|. Right: (Color online) Argument of 
the pQCD folding integral Eq. (|8j) on (y,y ma x) based on p-p 
fragmentation functions. 



xx, e(6rj, A?y) is the efficiency for detecting jet fragments 
in r\ acceptance br\ relative to 4-7T acceptance An, and 
do~ dijet /dy m ax is the dijet spectrum on parton rapidity. 

Figure [5] (right panel) shows the integrand of the fold- 
ing integral for p-p collisions. The p-p FFs are based 
on a parameterization of e + -e~ FFs with lower limit 
0.35 [27J, but measured FFs from p-p colli- 
require a higher cutoff at y m i n ~1.5 represented 
by the horizontal dotted line. The parton spectrum is 



ymin ^ 

sions 4 



where d 2 nf l / 'dy dr\ is the FD describing a spectrum hard 
component, D xx is the FF ensemble from collision system 



the solid curve in the left panel. The LHS of Eq. [8] is the 
solid curve in Fig. [5] (right panel). The comparison with 
p-p hard component (solid points) determines the 3 GeV 
parton spectrum cutoff to about 5%. The dotted curves 
represent ±10% variations in the cutoff energy. 

Figure [10] illustrates extension of the fragment distri- 
bution comparison to A-A collisions. The data are spec- 
trum hard components from five centralities of 200 GeV 
Au-Au collisions [HI . The bold dashed curve through p-p 
(points) and peripheral Au-Au (60-80%) data is the solid 
curve from Fig. [6] (right panel) corresponding to a par- 
ton spectrum cutoff of 3.0 ± 0.15 GeV and y m in ~ 1-6. 
Bold dotted curve H ee - me d is a reference correspond- 
ing to medium-modified e + -e~ FFs with y m i n reduced 
to 0.35 and parton spectrum cutoff reduced to 2.7 GeV. 
The dashed curve labeled 5 passing through central (0- 
12%) data corresponds to the same conditions but with 
ymin ~1.2. The spectrum hard component, even in more- 
central Au-Au collisions, reveals the parton spectrum 
cutoff unambiguously. The cutoff for central Au-Au col- 
lisions is reduced by about 10% compared to the 3 GeV 
observed for p-p collisions, resulting in a 50% increase 
in the dijet cross section because of the steep power-law 
spectrum [l4| . 
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FIG. 10: (Color online) Measured spectrum hard components 
Haa for five centralities from 200 GeV Au-Au collisions (bold 
curves of several colors) [l3l ] and 200 GeV NSD p-p colli- 
sions (points) [351 ] compared to calculated FDs for several 
conditions (vacuum, medium, e + -e~, p-p). The hatched re- 
gion at upper left estimates the uncertainty due to the soft- 
component Snn subtraction common to all centralities. 



C. Final-state particle, p t , E t production 

The relation between final-state hadron and p t / E t pro- 
duction and the initial partonic state falls between two 
extremes: 1) a two-component system with fragmenting 
free partons and parton spectrum with cutoff at 3 GeV 
(GLS) and 2) a single component of thcrmalized partons 
from a spectrum with cutoff near 1 GeV (SSM) transi- 
tioning to hadrons. The total particle, pt or E t produc- 
tion density is less informative than what is produced per 
participant nucleon pair in one unit of 77. How are N-N 
collisions in A-A different from isolated p-p collisions? 

1) GLS - In 200 GeV NSD p-p collisions the charged- 
hadron yield is dn c h/drj ~ 2.5, with hard-component con- 
tribution 0.02 derived from the n c h dependence [35]. The 
corresponding dP t /dri w 2.5 x 0.35 + 0.02 x 1 w 0.90 
GeV/c. 

In central 200 GeV Au-Au collisions the per- 
participant-pair hard-component yield increases relative 
to N-N collisions by a factor 54 = 6 (mean number of 
N-N collisions per projectile nucleon) x 3 (increase in di- 
jet multiplicity) x 1.5 (increase in minijet cross section) 
x 2 (increase in jet efficiency for many jets per colli- 
sion) [l4T ]. The charged multiplicity per participant pair 
is then (2/n par t) dn c hjdr\ = 2.5 + 1.1 = 3.6 (consistent 
with data from Au-Au central collisions), and the total-pt 
density is (2/n part ) dP t /drj = 2.5x0.35 + 1.1x0.45 = 1.38 
GeV/c since (pt) for the hard component is reduced from 



1 to ~ 0.45 GeV/c in central Au-Au collisions due to 
medium modification of FFs [3] • The estimates assume 
that the soft component remains unchanged with A-A 
centrality. There is no evidence to the contrary. 

2) SSM - An example of the single-component ap- 
proach estimates dnj et /dy = 750 in central Au-Au colli- 
sions assuming a parton spectrum cutoff at 1 GeV [20j. 
That parton (jet) density implies (2 / n part ) dPt / drj — 
2/3 x 750/191 x 1 GeV/c = 2.6 GeV/c assuming 2/3 of 
jet fragments are charged (pions) and there are 191 par- 
ticipant pairs in central Au-Au collisions. The pt density 
is twice what is observed in data, and the integrated dijet 
cross section is greater than the N-N total cross section. 

In another approach NLO first moment a(E t ) Po in one 
unit of rapidity for central Au-Au collisions is estimated 
to be 8 and 120 mb-GeV for parton spectrum cutoffs 
Po = 3 and 1 GeV respectively For T AA (b = 0) ~ 
34 mb" 1 [IH we obtain (2 / n par t) dE t / drj ~ 1.4 and 21 
GeV for po = 3 and 1 GeV respectively. A significant 
(tens of percent) reduction of those estimates would re- 
sult from excluding the Et < po region. LO estimates are 
~2 times smaller. Such estimates are typically compared 
with the entire final state, not a separate (semi)hard com- 
ponent explicitly identified with large-angle parton scat- 
tering. Data are consistent with {2 / n part ) dP t / drj ~ 0.5 
GeV/c for charged particles from observed minijets. 

Aside from the excess pt/E t compared with data 
the single-component approach ignores soft-component 
contributions to the final state in more-central A-A 
collisions without explaining what process causes the 
change. Peripheral Au-Au collisions follow the GLS ref- 
erence, including contributions from soft and hard com- 
ponents [131 ] . At what point does the dominant soft com- 
ponent disappear and its final-state contribution revert 
to scattered partons? Equivalently, what happens to the 
PDFs of projectile nucleons? 

IX. HYDRO AND SP SPECTRA 

A hydro description of nuclear collisions requires that 
radial flow be observed in p t spectra as a manifestation 
of large thermalized energy density and corresponding 
radial pressure gradients [3j. Evidence for radial flow in 
the final state is sought via blast wave (BW) fits to SP 
spectra [58l. [59|. 

A. p-p spectra 

Figure [Tl] illustrates attempts to describe p-p pt spec- 
tra with the blast-wave model (assuming hydro expan- 
sion in elementary collisions) [60]. In the left panel 
a simple procedure is illustrated to generate a BW 
model spectrum. A Maxwell-Boltzmann (MB) spec- 
trum 1/mt dn/dmt oc exp(— rrit/T) with slope parameter 
T = 0.145 GeV transformed to yt with mt = mo cosh(y t ) 
is boosted (blue shifted) by varying amounts Ay t ~ flt 
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according to Hubble expansion (boost proportional to 
radius), with maximum boost Ay t — 0.5. The result- 
ing BW spectrum corresponds to mean transverse speed 
(Pt) ~ 0.25 sometimes attributed to p-p collisions [60]. 
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FIG. 11: Left panel: Maxwell-Boltzmann (MB, exponential 
on m t ) distributions on transverse rapidity yt for a linear 
distribution of boosts Ay t representing the blast-wave (BW) 
model, with (Ayt) = 0.25 and T = 0.145 GeV. Right panel: 
Two-component model of the 200 GeV NSD p-p y t spectrum 
(dash-dotted curve), soft component Snn (dashed curve) and 
hard component Hnn (light solid curve). The bold solid curve 
is the BW model obtained by averaging the left panel on Ayt . 
The dotted curve is the MB limit of the BW for (Ay t ) -> 0. 

Figure [TTJ (right panel) compares the resulting BW 
model (bold solid curve) with the two-component rep- 
resentation of 200 GeV p-p data (dash-dotted curve) de- 
scribed in Sec. IV A} with soft component Snn (dashed 
curve) and hard component 0.02 Hq (light solid curve). 
BW fits are typically restricted to a "hydro" y t inter- 
val bounded above by p t ~ 2 GeV/c (y t ~ 3.3, ver- 
tical dashed line). BW parameter T is the same two- 
component value 0.145 GeV determined in (35[. (Ay t ) = 
0.25 w (Pt) is the mean "radial flow" value obtained from 
BW fits to p-p spectra (60|. 

Although the BW model appears to describe p-p spec- 
tra within a restricted yt/pt interval it cannot describe 
a more-extended y t interval because the curvature is 
much larger than typical data trends. There is also no 
absolute-yield constraint to BW fits. In contrast, the 
two-component spectrum model includes a GLS reference 
which specifies expected absolute magnitudes for p-p rich 
dependence and all A-A centralities in the event of N-N 
linear superposition and binary-collision scaling. 

This comparison indicates that for p-p spectra the BW 
model accommodates the Levy form of soft component 
SNN{yt) which describes the zero-density limit of p-p col- 
lisions where radial flow would be least likely. 



B. A-A spectra 

An early application of the BW model to S-S spec- 
tra is discussed in Sec. IV Bl The S-S result is iden- 
tical to the 200 GeV p-p case because the spectrum is 
well-described by Snn for p-p collisions. S-S collisions 



at ^/snn = 19.4 GeV appear to be well-described by a 
200 GeV GLS reference with negligible hard component. 
The BW model has been applied more recently to RHIC 
Au-Au spectra to infer systematic variation of transverse 
speed (f3t) and kinetic decoupling temperature Tki n [H|. 
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FIG. 12: Left panel: Similar to Fig. [TT] (right panel) but for 
200 GeV central Au-Au collisions with v = 6. The BW pa- 
rameters for this case are {Ayt) = 0.6 and T = 0.10 GeV. 
Right panel: {(3 t ) data for 62 and 200 GeV Au-Au colli- 
sions [611 ] vs centrality measure v. Construction of lines and 
curves is described in the text. The hatched region locates a 
sharp transition in jet correlation characteristics in 200 GeV 
Au-Au collisions. To the left of the transition jet correlations 
and yields follow binary-collision scaling of p-p collisions. To 
the right jet yields increase dramatically [161 ]. 

Figure [T2] (left panel) shows the two-component ref- 
erence (y = 6, dash-dotted curve) for central Au-Au 
collisions at 200 GeV and corresponding fixed soft com- 
ponent Snn- Although spectrum data deviate substan- 
tially from the two-component reference at larger y t com- 
parisons with hydro near p t ~ 1.5 GeV/c (yt ~ 3) are 
meaningful. The bold solid curve is the BW model with 
{(3 t ) = 0.6 and T^ in = 0.10 GeV, values corresponding to 
typical evolution of fitted BW parameters with HI colli- 
sion centrality as shown in Fig. [T2l fright panel). Within 
the restricted yt interval conventionally assigned to hy- 
dro (left of the dashed line) the BW model again seems 
to describe the data well, but fails elsewhere. The BW is 
determined for central Au-Au collisions primarily by the 
spectrum hard component (light solid curve) attributed 
to jets. 

Figure [T2] (right panel) shows published (f3t) (points) 
from BW fits to identified-particle spectra for p-p colli- 
sions at 200 GeV [6fJ and Au-Au collisions at ^/snn = 62 
and 200 GeV [6l[ plotted on mean participant pathlength 
v [26j . The solid curve for 200 GeV is constructed as fol- 
lows. The dashed line is a fit to the right-most six points. 
The dash-dotted line passing through the p-p point cx v 
represents a Glauber linear superposition (GLS) refer- 
ence for minijet production. The solid curve then inter- 
polates the lines. The hatched region denotes a sharp 
transition observed in 200 GeV minijet characteristics 
from correlation data flU ] and in the Au-Au spectrum 
hard component for pions and protons [l3| , visible as the 
nonuniform trend with gap at large yt in Fig. [2] (right 
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panel) . The 62 GeV dotted curve is similarly constructed 
and also corresponds to jet correlations. 

The BW model appears to describe spectra in the re- 
stricted interval p t <E [0.5,2] GeV/c (y t <E [2,3.3]). Com- 
parison with the two-component model shows that evolu- 
tion of BW parameters «/3 t ),T fei „) from (0.25,0.15 GeV) 
to (0.6,0.09 GeV) corresponds to at least six-fold increase 
of the hard component with centrality (cf. Fig. [T2")) . BW 
parameter variations are seductive because they suggests 
isentropic expansion of a thermalized hadron fluid, with 
increasing separation between "chemical" and "kinetic" 
decoupling temperatures Q- However, the BW competes 
with QCD processes already identified in elementary col- 
lisions and expected to follow binary-collision scaling in 
A- A collisions. Contrast the BW description of A- A spec- 
tra with details of hard-component evolution revealed by 
taa in Fig. [2] (right panel). The "temperature" parame- 
ter of soft component Snn (equivalent to Tkinetic) shows 
no change with A-A centrality. 

This reinterpretation of BW trends in terms of frag- 
mentation systematics (jet correlations) illustrates the 
importance of the two-component fragmentation refer- 
ence for nuclear collisions. As a rule, QCD mechanisms 
should be given first option to describe data before hydro 
conjectures are imposed. 



C. Identified-proton spectra challenge radial flow 

Figure [T3] shows proton spectra from Au-Au colli- 
sions for five centralities at ^Jsnn = 200 GeV (five 
solid curves) [l3j]. The 3D density has the form p p — 
{l/2iryt) d 2 n p /dy t dy z . Soft- and hard-component refer- 
ence functions Snn and Hnn are defined as asymptotic 
limits of Au-Au proton spectrum v dependence consis- 
tent with a similar analysis of p-p n c h dependence [35| . 
The dashed curves represent GLS reference Snn+v Hnn 
for limiting cases v — 1 and 6. "Jet quenching" suppres- 
sion relative to the reference is apparent above y t ~ 4.5 
(p t ~ 6 GeV/c). The excess proton yield near y t = 3.7 
(p t ~ 2.5 GeV/c) corresponds to the baryon/meson "puz- 
zle" at RHIC. Because of their greater mass protons 
(baryons) should be more sensitive to a boosted hadron 
source (radial flow). 

Arrows and labels at the bottom of the figure indicate 
conventional division of pt spectra into hydro, ReCo and 
pQCD intervals. The corresponding p t values are indi- 
cated at the upper margin. As noted, hydro phenomena 
are expected to dominate for p t < 2 GeV/c (y t < 3.3). 
However, proton spectra below y t ~ 3 (p t ~ 1.5 GeV/c) 
closely follow the GLS reference. Near y t ~ 2.7 hard 
component Haa dominates the proton spectrum and fol- 
lows binary- collision scaling to the error limits of data. 
That trend is not expected if radial flow drives Au-Au 
spectra as assumed in the BW model. 

If radial flow played a significant role yt spectra should 
be boosted increasingly to the right (positive Ay to) with 
increasing centrality (and radial pressure), which trend 




y t (n) 

FIG. 13: (Color online) Proton yt spectra for five Au-Au 
centralities (solid curves) The dashed curves are two- 

component model functions for N-N {y = 1) and b = (y = 6) 
Au-Au collisions. The dash-dotted curve is the N-N hard- 
component reference function. The dotted curve labeled Snn 
is the soft-component model for protons from Au-Au colli- 
sions. 



is not observed. In this direct confrontation with hy- 
dro, proton data reveal that the majority of final-state 
hadrons does not manifest radial flow in Au-Au collisions 
at RHIC. Instead, a large fraction of protons emerges 
from parton fragmentation. Corresponding pion spectra 
provide similar conclusions [l3j]. 



X. HYDRO AND AZIMUTH CORRELATIONS 

Just as evidence of radial flow is required for a hy- 
dro interpretation of heavy ion collisions, measurement 
of elliptic flow in non-central collisions (reflecting the ec- 
centric interaction region) is also essential for hydro. The 
conventional elliptic flow measure is V2 — (cos[2((j>— \Ev)]) 
measured relative to the reaction plane defined by the 
A-A impact parameter and collision axis [9]. Because 
the reaction plane is not observable vi analysis relies on 
variants of multi-particle azimuth correlations denoted 
by v 2 {method} [13, SI- 

In a hydro description the transverse-rapidity yt dis- 
tribution of a thermalized expanding medium would de- 
scribe a common boosted source for all hadron species. 
Radial boost distribution Ay t (r, <f>) is then equivalent to 
flow field f3t(r,<f>) [fit = tanh(Ay t )]. If the final-state 
hadron distribution is described by density p(yt , 4>, b) pt- 
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differential 1*2 (pt) is given by 



V2(Pt,b) 



J d(pp(p t ,(j>, b) cos(20) 
J d(f)p{jp t ,(f),b) 



(9) 



with 4> defined relative to the reaction plane. p t -integral 
V2(b) is obtained by integrating numerator and denomi- 
nator over p t . V2(b) is compared to eccentricity e(b) of the 
initial A-A configuration space as a test of hydro models. 
Although the V2 concept appears simple, implementation 
and interpretation of «2 analysis is not. 

Strong minijet contributions to azimuth correlations 
(called "nonffow" ) dominate V2 data obtained with con- 
ventional methods over a substantial part of the A-A 
centrality interval [25| . making interpretation of V2 data 
within the hydro context problematic. The model- 
neutral term "azimuth quadrupole" can be applied to w 2 
as a measure of azimuth correlations. Quadrupole corre- 
lations are indeed observed in RHIC data, but although 
v 2 is conventionally interpreted as elliptic flow [24[ phys- 
ical interpretation of the quadrupole is an open ques- 
tion [lj| 0, (25[. Absence of measurable radial flow in 
SP spectra contradicts the concept of elliptic flow as a 
modulation of radial flow in non-central collisions. 



A. pt-integral V2 

The centrality and energy systematics of ratio i^/e 
have been interpreted to indicate monotonic increase 
with a (Knudsen) number of secondary collisions (low- 
density limit or LDL model) suggesting approach to ther- 
malization [62| . Saturation of V2I £ could indicate a hydro 
limit. Claims for achievement of the hydro limit [631 ] and 
so-called constituent-quark scaling (64[ Hfl have been in- 
terpreted as evidence for a thermalized QGP at RHIC. 

Fig. [T3] (left panel) shows per-particle azimuth 
quadrupole measure Ap[2]/^/p ro f on centrality measure 
v [26]. The quadrupole amplitude can be obtained from 
fits to 2D angular autocorrelations 0, H3, [Hj], and is 
related to V2 by 



Ap[2] 



dr] d(f> 



[Mm? 



(10) 



defining v 2 {2D} (TaHl^l- The definition is compatible 
with minijet angular correlation measurements and the 
Glauber linear superposition (GLS) reference [l6t | . 

Fig. [TJ] shows data (points) for conventional elliptic- 
flow measures t>2{2} and ?^_{4} (two- and four-particle 
correlations respectively) [66j. The solid curve (trans- 
formed from the straight line in the right panel) is defined 
by 



= 0.0045 i?(Vs^)C t n bm . (11) 



R = 111(7^/13.5 GeV)/ ln(200 GeV/13.5 GeV) is 
consistent with the energy dependence of pt correlations 
attributed to minijets [181 ]. The dotted curve passing 
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FIG. 14: Left panel: pj-integrated V2(b) data (solid points, 
open triangles) for two analysis methods from 200 GeV Au- 
Au collisions vs centrality measure v [§6j]. The open squares 
are v 2 {EP} data from 17 GeV Pb-Pb collisions [63] ■ The bold 
solid curve is derived from the line from the right panel. The 
dotted curve includes jet correlations (same-side jet peak). 
The light solid and dashed curves apply to the 17-GeV data. 
Right panel: Data and curves from the left panel divided by 
Optical- The solid line is defined by Eq. (|11[) . The dash-dotted 



curve indicates the trend if t Monte 



were used instead. 



through f2{2} data is obtained by adding to Eq. (fTTj) 
"nonflow" term g2/2ir = 0.004 u 5 representing the m = 
2 Fourier component of the measured same-side jet peak 
on azimuth [l| (cf. Sec. EE]) , 

For 17-GeV Pb-Pb v 2 {EP} (event-plane) data (open 
squares) [63] the same fractional relation between "flow" 
(thin solid curve) and "nonflow" (dashed curve, with 52 
added) is observed. Thus, Eq. (jll[) with nonflow param- 
eter 32 obtained from (mini)jet measurements describes 
all Au-Au (or Pb-Pb) V2 data above ^/snn — 13.5 GeV. 

Fig. [14] (right panel) shows the same data in a different 
format. The optical eccentricity e op t leads to a simple, 
universal linear relation. Eq. (|lip defines the solid line. 
An alternative Monte Carlo eccentricity would produce 
the dash-dotted curve. Arguments supporting the optical 
eccentricity are given in (25|. Deviations of the dotted 
and dashed curves represent the effect of jets on V2{EP} 
and «2{2} measurements. 

A combination of V2 data, minijet correlations and the 
optical-Glauber description of collision eccentricity re- 
veals that over a large range of energies and centralities 
V2 data are described accurately and completely by the 
initial state (6, ^snn) [HL [25[. There is no apparent 
dependence on subsequent collision evolution, equation 
of state or degree of thermalization. There is no cor- 
respondence with the recently-observed sharp transition 
in minijet angular correlations [l6| , calling into question 
any relation among the azimuth quadrupole, conjectured 
bulk medium and hydrodynamics. 



B. pt-differential V2 

Fig. [T5l fleftpanel) shows V2(pt) data for three hadron 
species [32], [33[. The mass ordering below 2 GeV/c is 
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interpreted to imply a hydro phenomenon. Baryon vs 
meson trends at larger p t are interpreted to conclude 
that constituent quarks coalesce to form hadrons which 
exhibit elliptic flow, implying a thcrmalized flowing par- 
tonic medium prior to decoupling (cf. Sec. IX C"]) . Dotted 
theory curves A and B for pions are from [68j| and [6^ | re- 
spectively. The curves passing through data are from (l9j 
(and cf. Fig. [TSJ). 




p t (lab) (GeV/c) y,(7i,K,p) 



FIG. 15: Left panel: t>2(pt)data for three hadron species 
plotted in the usual format [3^. |33||. Right panel: The same 
V2(pt) data divided by p t in the lab frame suggest universal- 
ity on proper transverse rapidity for each hadron species — in 
particular the correspondence of data near yt = 1. Dotted 
curves A and B in each panel are viscous hydro predictions 
from [68|] and respectively. The three curves through data 
are derived from this analysis. 

Fig. [15] (right panel) shows the same data plotted as 
ratio Vijpt vs yt(^,K,p) denoting transverse rapidity (a 
velocity measure) computed with the correct mass for 
each hadron species. Motivation for the v^jvt ratio is 
described in Sec. lIVCl The three hadron species follow a 
common trend at smaller y t expected for hadrons emitted 
from a common boosted source. The mean boost Ay t ~ 
0.6 is approximately the source radial speed, since (3 t = 
tanh(Ay t ). Although a hydro mechanism might produce 
such a boost, these data to not require a hydro boost 
mechanism. 

Theory curves A and B are zero- viscosity limits of two 
viscous- hydro calculations [68|, [69[ . The differences are 
most notable in the plot on the right. At small y t the 
theory boost distributions, a central element of the hydro 
model, are very different from data. The hydro calcula- 
tions are consistent with radial Hubble expansion of a 
fluid, whereas the data are consistent with an expanding 
cylindrical shell [l9[ . 

Fig. [TS] shows quadrupole spectra (points) recon- 
structed as described in Eq. ([3} and [19( from the 
V2(pt) data plotted in Fig. [15] The additional factor 
{2/n part ) p(yt) relative to Fig. [T5] (right panel) is ob- 
tained from parametrizations of single-particle identificd- 
hadron spectra in (l3j |. The SP spectrum factor elim- 
inates the extraneous denominator in Eq. ([9} to reveal 
the quadrupole spectrum implicit in the numerator. All 
aspects of minimum-bias v^ipt) structure are simply and 
accurately described. 



10 



10 







i i i | i i i i | i i i i | i i i 








^NNp 




minimum-bias 




200GeVAu-Au 
















A \ > 


v - 3.5 

\ , 


; 


.\\ \ V, K 




; 




boosted \ 




A 2 - 0.005 \ \ \ 
T, -O.lGeV *\ , 


hadron \ 
source \ 




Ay 10 = 0.6 \ \ 
"\ ■ 


, ,\ , i , , , , i , , , , \ , , A" 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



y,(7i,K,p) 

FIG. 16: Data from Fig. [15] transformed to proper y t for each 
hadron species [IsJ. The dotted curves are soft components 
from respective single-particle spectra for comparison. The 
prominent feature is the common edge at y t ~ 0.6, imply- 
ing that the quadrupole components for three hadron species 
originate from a common boosted source. Hadron abundances 
and quadrupole spectrum shapes are similar to the single- 
particle spectrum soft components. 



The quadrupole spectra are shifted by common 
monopole boost Ay t0 ~ 0.6 (cf. the left edges). The 
spectrum shapes (solid curves) have the same Levy dis- 
tribution form on nit that describes the SP soft com- 
ponents (Snnx, dotted curves), with the same relative 
hadron abundances but substantially different parame- 
ters. The "temperatures" of the quadrupole spectra are 
systematically reduced by at least 30%, e.g. to T2 ~ 0.10 
GeV from the SP soft-component T = 0.145 GeV for 
pions. The spectrum widths on y t of quadrupole spectra 
vary with ratio T 2 /mo just as f° r SP spectra. The solid 
curves, transformed to other plotting formats in Figs. 1151 
and 1 171 describe quadrupole data within errors. 

Quadrupole spectra are reduced in amplitude relative 
to SP spectra by a common factor A2 ~ 0.005. The 
reduction factor includes the product of the true mo- 
mentum eccentricity represented by quadrupole boost 
Ayt2 and the quadrupole spectrum absolute yield. Ar- 
guments in [19( suggest that Ayt2 ~ Ayto and the actual 
quadrupole yield is a small fraction (< 5%) of the total 
multiplicity. 

The information derived from V2(pt) data is thus a 
boost distribution (narrow relative to the mean, inconsis- 
tent with Hubble expansion), a quadrupole source "tem- 
perature," a relative abundance, and distinction between 
SP spectrum soft component and quadrupole component. 
The combination of properties is inconsistent with a ther- 



16 



malized fluid medium as the source for most hadrons from 
a collision. Hadronization appears to follow the same pro- 
cess for soft and quadrupole spectrum components, but 
the hadron sources are distinct. 

The sequence of three panels in Figs. [15] and [TH] re- 
veals that v 2 (pt) as defined exaggerates the role of data 
at larger p t due to the presence of the single-particle spec- 
trum in its denominator and a kinematic factor p t in its 
numerator [l9|. The most important v 2 data from the 
standpoint of the hydro model are at smaller pt, espe- 
cially for low-mass hadrons, as is apparent from Fig. 1161 
Those details are suppressed by the conventional plotting 
format of Fig. [15] (left panel). 

C. Constituent-quark scaling 

The same v 2 (pt) data can be replotted in a format used 
to demonstrate so-called "constituent quark scaling" [64|, 
Hlj]. Fig. [T7J (left panel) shows t>2 (p*) plotted vs transverse 
mass nit — y/ Pt + JWq . Offsets due to the common source 
boost at smaller p t are reduced on wit, and there may be 
vertical segregation into mesons and baryons at larger 




m t -m (GeV/c 2 ) 2(m t -m )/n q (GeV/c 2 ) 



FIG. 17: Left panel: V2(pt) data from Fig. [15] are plotted vs 
kinetic energy mt — mo. The intercept spacing at small pt 
is reduced by a factor 3x. Right panel: The left panel, but 
with constituent quark scaling in the form 2/n q so that meson 
trends remain unchanged. The shift for baryons is indicated 
by the arrows. 

Fig. [13 (right panel) is said to demonstrate constituent- 
quark scaling. The number of valence quarks n q for 
hadrons (2 or 3) is applied to data. Scaling implies that 
data for all hadrons should lie on a common curve, which 
seems to be approximately the case below 2 GeV/c 2 . Ad- 
ditional factors 2 in this plot are included so that only 
baryon data shift from left to right panel. The curves 
through data from pjl and Fig. [TH] are transformed ap- 
propriately. The baryon data seem closer to meson data 
in the right panel within a limited mt interval. 

Do those results support constituent-quark scaling? 
The boosted intercept points on p t in Fig. [15] (left panel) 
fall at p t0 = m sinh(Ay t0 ) ~ m Ay m , with Ay m = 
0.6. The intercepts scale with hadron mass, and the 
common boost can thus be read off the plot directly. 



The corresponding intercepts in Fig. [17] (left panel) are 
(m m - m ) = mo [cosh(Ay t0 ) - 1] ~ m (Ay m ) 2 /2. Sen- 
sitivity to the common boost has been reduced by a fac- 
tor Ay t o/2 ~ 1/3, minimizing the hydro phenomenon in 
question relative to data errors. The curves from p^ |. 
which describe data well, diverge from one another above 
2 GeV/c 2 and still retain significant spacings at small mt 
due to the common source boost. 

In Fig. [17] (right panel) additional factors n q /2 have 
proportions 1:1:1.5 for the three hadron species, whereas 
v^ipt) hi all its aspects depends on hadron masses in ap- 
proximate proportions 1:3.5:7. The boost manifestation 
on momentum, the quadrupole spectrum shapes (widths) 
and the fractional abundances are all determined by 
hadron masses, most evidently in Fig. 1161 Since hadron 
masses don't scale as the number of valence quarks the 
constituent-quark scaling hypothesis is falsified by V2(pt) 
data, despite the appearance of Fig. [T7J (right panel). 

Fig. |T5] (right panel) provides the strongest chal- 
lenge to constituent-quark scaling. The measured 
quantities appear as ratios — v^/pt and ratios within 
ytipt/fno^mt/mo) — for which constituent-quark num- 
bers should cancel, and where quark coalescence in a 
common boosted frame should be most naturally de- 
scribed. There is no scaling. The large differences 
between hadron species are simply explained by mass- 
dependent effects in the quadrupole spectra of Fig. [TB] 

XI. DISCUSSION 

Results from hydro-motivated analysis methods are 
conventionally interpreted to favor hydro mechanisms 
and a dense partonic medium. Is formation of an sQGP 
with anomalously small viscosity required by such results? 
Do other results contradict such conclusions? 



A. Minijets, initial conditions and thermalization 

Thermalized minijets are expected to provide the large 
energy densities required to drive hydro expansion in 
more-central RHIC collisions 0, H, 0, [48[. Evidence 
for partonic hydro expansion would in turn confirm for- 
mation of a thermalized QCD medium or QGP arising 
from parton (minijet) multiple scattering. High-p t jet 
tomography would map the resulting "opaque core" and 
determine its properties [7(| • The opaque core 0] is 
said to reduce the sensitivity of ratio measure Raa due 
to the predominance of "surface" jets emitted from a thin 
"corona" layer [7lj . At most 1% of hadrons with up to 14 
GeV/c momentum from "quenched jets" would survive 
according to 0. Those expectations are strongly con- 
tradicted by some data, as demonstrated in this analysis. 

The most important challenge to hydro at RHIC is 
the survival of essentially all minijets to the final state 
in central 200 GeV Au-Au collisions EE EE El 0- 
Minijet analysis of spectra and correlations reveals that 
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even central Au-Au collisions are nearly transparent to all 
scattered partons (not just "surface" partons). Back-to- 
back jet correlations are not diminished. There is no re- 
duction of correlations by parton or hadron rescattering. 
Hadron fragments down to 0.1 GeV/c retain their angu- 
lar correlations relative to parent partons. However, par- 
ton fragmentation is strongly modified in more-central 
collisions [141 ] . 

B. Final-state particle and p± production 

Some arguments for parton thermalization have been 
based on a conjectured saturation-scale parton spectrum 
cutoff near 1 GeV which would increase the scattered- 
parton yield and phase-space density 10-30 fold (com- 
pared to a 3-GeV cutoff). Such large densities could 
insure thermalization by parton multiple scattering, but 
they seem to be inconsistent with the observed final state. 

A quantitative relation has been established 
among parton cross sections, FFs and hadron spec- 
tra/correlations [3]. Comparison of pQCD elements 
to p-p hadron spectra establishes a parton spectrum 
cutoff near 3 GeV. Peripheral Au-Au collisions follow a 
two-component GLS reference up to a transition point 
on centrality. Beyond that point spectra and correlations 
are still described by parton FFs and a parton spectrum 
terminating near 3 GeV, although there are significant 
changes to fragmentation. The relation between soft and 
hard spectrum components follows pQCD expectations. 

In central Au-Au collisions the soft component (projec- 
tile nucleon fragmentation) represents at least two-thirds 
of hadron production. The other third is accounted for by 
additional parton scattering and fragmentation compared 
to the GLS reference, with modified FFs. pt/E t produc- 
tion is similarly well described. Essentially all semihard 
scattered partons appear as correlated hadrons in the fi- 
nal state, contradicting claims of parton thermalization 
and opaque medium. 



C. Radial flow 

Radial flow should be the primary manifestation of 
thermalized initial-state energy densities and pressure 
gradients. Radial flow inferred from hydro-inspired blast- 
wave fits to hadron spectra seems to confirm required 
large energy densities. But conventional BW models as- 
sume that any deviation from a Maxwell-Boltzmann dis- 
tribution below 2-3 GeV/c is caused by radial flow. 

The BW spectrum model must compete with the two- 
component model and contributions from minimum-bias 
parton fragmentation. The two-component model cou- 
pled with Glauber linear superposition describe data fea- 
tures associated with known QCD processes which should 
a priori appear in A-A collisions: projectile nucleon frag- 
mentation and large-angle parton scattering and frag- 
mentation consistent with elementary collisions. The 



BW model approximates spectra with limited accuracy 
over a limited p t interval and fails qualitatively elsewhere. 
In effect, parton fragmentation processes are misidenti- 
fied by BW fits and injected into the hydro context. 



D. Elliptic flow 

V2 interpreted as elliptic flow is seen as the most con- 
vincing support for hydrodynamics at RHIC. The v 2 def- 
inition assumes a monolithic flowing medium common to 
most final-state hadrons. But quadrupole spectra ex- 
tracted from published ^(p*) data are quite different 
from corresponding SP spectra. Only a small fraction 
of the hadronic final state may actually "carry" the az- 
imuth quadrupole [19J. 

Sharp transitions in minijet characteristics at specific 
A-A centralities [l6| could be interpreted to indicate 
medium modification of parton fragmentation. v 2 data 
show no evidence of such transitions, in conflict with their 
interpretation in terms of elliptic flow of a homogeneous 
thermalized medium. 

The concept that elliptic flow results from rescatter- 
ing (of partons or hadrons) in a medium is contradicted 
by minijet correlation data. Such rescattering would de- 
couple partons from their initial-state scattering partners 
(especially at 3 GeV) and destroy hadronic correlations 
resulting from parton fragmentation. Strong back-to- 
back parton correlations are inferred from strong hadron 
jet correlations, even in central Au-Au collisions. 



XII. SUMMARY 

Differential analysis of single-particle spectra and two- 
particle correlations from RHIC Au-Au collisions reveals 
that copious parton fragment yields extend down to small 
transverse momentum even in central collisions. Accord- 
ing to pQCD and elementary collisions at least half of 
all parton fragments in nuclear collisions should appear 
below 1 GeV/c. Measured hard components of spectra 
and correlations at RHIC are consistent with that expec- 
tation. Hard components are identified as fragment dis- 
tributions from minimum-bias parton fragmentation or 
"minijets" quantitatively described by pQCD. Observ- 
able minijets play a central role in RHIC collisions. 

In the hydro context single-particle spectra are divided 
into three regions, with the interval pt < 2 GeV/c as- 
signed to hydro phenomena. But that interval contains 
the peak of the spectrum hard component and most par- 
ton fragments. The blast-wave model fitted to SP spec- 
tra returns an estimate of radial flow. But the spectrum 
structure responsible is just the hard component associ- 
ated with minijets. Spectrum ratio Raa used to estimate 
jet suppression at larger p t strongly suppresses spectrum 
structure below p t ~ 4 GeV/c, including most of the jet 
fragment yield in the spectrum hard component. 
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Conventional V2 analysis, especially in more-peripheral 
and more-central A-A collisions, misidentifies the rn = 2 
Fourier component of jet azimuth structure ("nonflow") 
as elliptic flow and can greatly overestimate vi in such 
cases. Ratio v^ipi) contains the single-particle spectrum 
in its denominator, and is thus strongly affected by the 
spectrum hard component (minijets) at larger p t , in ad- 
dition to nonflow contributions from jet azimuth correla- 
tions in its numerator. 

So-called "triggered" analysis of jet azimuth correla- 
tions suppresses true jet yields because of incorrect es- 
timation of the combinatoric-background offset (ZYAM 
estimate) and v-i component. The conventional ZYAM 
procedure underestimates jet yields by as much as a fac- 
tor 10, and the away-side jet is distorted by the vi over- 
subtraction (minimum at n) so as to suggest the presence 
of Mach shocks. The combination leads to inference of 
an "opaque medium" and parton thermalization. 

The numerator of V2(j>t) contains quadrupole spectrum 
P2(yt) which could represent a hydro phenomenon if that 
were relevant. However, the quadrupole component is 
insensitive to dramatic changes in jet properties with in- 



creasing centrality that could be attributed to a common 
medium. The quadrupole component seems to be carried 
by an isolated and small part of the final-state system. 

In conclusion, the combination of two fragmentation 
components (projectile nucleons and semihard-scattered 
partons) plus an isolated third (quadrupole) component 
accurately describes all RHIC spectrum and correlation 
data over a large pt interval. No evidence for radial 
flow is observed. Jet yields increase according to binary- 
collision scaling in more-peripheral collisions, and even 
more rapidly for more-central Au-Au collisions. Frag- 
mentation is strongly modified in the latter case. There 
is no evidence for an opaque medium. 

As a rule, above ^/snn ~ 15 GeV QCD should be 
given first opportunity to describe nuclear collisions. Hy- 
dro should not be invoked until QCD has been shown 
to fail. A hydro description of selected aspects of data 
within restricted kinematic intervals does not establish 
the necessity of the model. In some cases hydro is con- 
tradicted. The relevance of hydrodynamics to RHIC col- 
lisions can therefore be questioned. 
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